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Using an isospin- and momentum-dependent modified Gogny (MDI) interaction, the 
Skyrme-Hartree-Fock (SHF) approach, and a phenomenological modified Skyrmc-like 
(MSL) model, we have studied the incomprcssibility X sa t(<5) of isospin asymmetric nu- 
clear matter at its saturation density. Our results show that in the expansion of A" sa t (<5) 
in powers of isospin asymmetry <5, i.e., X sa t(<5) = Kq + -fST S at,2<5 2 + -fsat,4<5 4 +0(<5 6 ), the 
magnitude of the 4th-order K aat ,4, parameter is generally small. The 2nd-order K sa , t ,2 
parameter thus essentially characterizes the isospin dependence of the incompressibil- 
ity of asymmetric nuclear matter at saturation density. Furthermore, the K aa x,2 can be 
expressed as Ksat,2 = Ksym — 6L — m terms of the slope parameter L and the 

curvature parameter K Byln of the symmetry energy and the third-order derivative pa- 
rameter Jo of the energy of symmetric nuclear matter at saturation density, and we find 
the higher order Jo contribution to -K" S at,2 generally cannot be neglected. Also, we have 
found a linear correlation between K Byln and L as well as between Jq/Kq and Kq. Us- 
ing these correlations together with the empirical constraints on Kq and L, the nuclear 
symmetry energy -E S ym(po) a t normal nuclear density, and the nucleon effective mass, 
we have obtained an estimated value of K Ba ,t,2 = —370 ± 120 MeV for the 2nd-ordcr 
parameter in the isospin asymmetry expansion of the incomprcssibility of asymmetric 
nuclear matter at its saturation density. 



1. Introduction 

With the establishment or construction of many radioactive beam facilities around 
the world, such as the Cooling Storage Ring (CSR) facility at HIRFL in China 
the Radioactive Ion Beam (RIB) Factory at RIKEN in Japan ^1, the FAIR/GSI in 
Germany^, SPIRAL2/GANIL in Fr ance^, and the Facility for Rare Isotope Beams 



1 
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(FRIB) in USA ^3 it is possible in terrestrial laboratories to explore the equation of 
state (EOS) of an isospin asymmetric nuclear matter under the extreme condition 
of large isospin asymmetry and thus to determine the density dependence of the 
nuclear symmetry energy. This knowledge is important for understanding not only 
the structure of radioactive nuclei, the reaction dynamics induced by rare isotopes, 
and the liquid-gas phase tra nsition in asymmetric nu clear matter, but also many 
critical issues in astrophysicsl 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 

For symmetric nuclear matter with equal fractions of neutrons and protons, 
its EOS is relatively well-determined. In particular, the incompressibility of sym- 
metric nuclear matter at its saturation density po has been determined to be 
240 ± 20 MeV from ana lyses of the nuclear giant monopole resonances (GMR) 
| 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | and itg EQS at densities of 2p < p < 5p has also been 

cons traine d by measurements of collective flows ^ and subthreshold kaon produc- 
tion^!^ in relativistic nucleus-nucleus collisions. On the other hand, the EOS of 
asymmetric nuclear matter and the density dependence of the nuclear symmetry 
energy, is largely unknown except at po where the nuclear symmetry energy has 
been determ ined to be around 30 MeV from the empirical liquid-drop mass for- 
mula, The values of the symmetry ene rgy at other de nsities, particularly at 

supra-saturation densities, are poorly known [6| 7 | 15 | 29 | 30 | 31 | 

Heavy-ion collisions, especially those induced by neutron-rich nuclei, provide a 
unique tool to investigate the EOS of asymmetric nuclear matter, especially the den- 
sity dependence of the nuclear symmetry energy. During the last decade, significant 
progress has been made both experimentally and theoretically on constraining the 
behavior of the symmet ry energ y at subsaturation density using heavy-ion reactions 
|32 | 33|34 | 35 | 36|37 | 38 | 3 l (g ee R e f.[H]f or the most recent review). More recently, there 
has also been an attempt SE1 to extract the symmetry energy at supersaturation 
densities from the FOPI data on the tt~ / tt + ratio in central heavy- ion collisions at 
sis/gsiHU. Information on the density dependence of the nuclear symmetry energy 
has also been obt aine d from the structure of finite nucl ei and their excitations, s uch 
as the mass data^Z] neutron skin in heavy nuclei ^2 giant dipole resonances 231, 
pygmy dipole resonance and so on. These studies have significantly improved 
our understanding of the EOS of isospin asymmetric nuclear matter. 

The incompressibility of asymmetric nuclear matter at its saturation density 
is a basic quantity to characterize its EOS. In principle, this informati on can be 
extracted experimentally by measuring the GMR in neutron-rich nuclei I 46 | 47 | j n 
the present talk, we report our recent work o n the i ncompressibility of isospin 
asymmetric nuclear matter at saturation density I 48 | 49 | 



2. Incompressibility of isospin asymmetric nuclear matter 

The EOS of isospin asymmetric nuclear matter, given by its binding energy per 
nucleoli, can be generally expressed as a power scries in the isospin asymmetry 
5 = (p n —Pp)/p, where p = p n + Pp is the baryon density with p n and Pp denoting the 
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neutron and proton densities, respectively. To the 4th-ordcr in isospin asymmetry, 
it is written as 

E(j>, 6) = E {p) + E sym (p)6 2 + E symA {p)S 4 + 0(5 6 ), (1) 

where Eq(p) = E{p, 8 = 0) is the binding energy per nucleon of symmetric nuclear 
matter, and 

l d 2 E(p,S) l &E{p,6) 

ti sym {p) - ^ -7^ |<5=0, &symA{P) = ggi l<5=0- 

In the above, E sym (p) is the so-called nuclear symmetry energy and E symA (p) is 
the 4th-ordcr nuclear symmetry energy. 

Around the normal density po, the E$(p) can be expanded, e.g., up to 4th-ordcr 
in density, as 

E (p) = E Q ( P0 ) + L 0X + |^x 2 + | X 3 + | X 4 + 0(x 5 ), 

where x is a dimcnsionlcss variable characterizing the deviations of the density 
from the saturation density p of symmetric nuclear matter, and it is conventionally 
defined as x — B j^ L - The Eq(po) is the binding energy per nucleon in symmetric 
nuclear matter at the saturation density po and the coefficients of other terms are 
defined by 

dE {p) _ 2 d 2 E (p) 

— dp — >p=p°> 0= — ^2 — \p=poi \ 6 ) 



r _g 7n 3 ^o(p) , r _ o ,A d A E {p) 

J n - tlpo ^3 |p=p > J Q — ° L Po \P=Po- 



E sym (p) = E sym (p ) + L X + ^x 2 + ~^X 3 + -^X 4 + 0(x% (5) 



Obviously, we have Lq = according to the definition of the saturation density po 
of symmetric nuclear matter. The coefficient Kq is the so-called incompressibility 
coefficient of symmetric nuclear matter, and it characterizes the curvature of Eo(p) 
at po. The coefficients Jo and Iq correspond, respectively, to 3rd-order and 4th-order 
incompressibility coefficients of symmetric nuclear matter. 

Around po, the E sym (p) can be similarly expanded up to 4th-order in x as 

^sym 2 , ^sym 3 , 4ym 4 
— - — Y H — Y H — 

2! 3! X 4! 

where L, K aym , J sym and I sym are the slope parameter, curvature parameter, 3rd- 
order coefficient, and 4th-order coefficient of the nuclear symmetry energy at po- 
We can also expand the 4th-order nuclear symmetry energy E syuiA (p) around po 
up to 4th-order in x as 

E symA (p) = E symA (Po) + L symAX + + ^gpV + ^X 4 + 0(x% (6) 

with L synii 4, i^ symi 4, J S ym, 4 and 4ym,4 being the slope parameter, curvature pa- 
rameter, 3rd-order coefficient, and 4th-order coefficient of the 4th-order nuclear 
symmetry energy E symA (p) at p . 
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Conventionally, the incompressibility of asymmetric nuclear matter is defined 
at its saturation densi ty p sa ,t, where P(p, S) = 0, thus also called the isobaric incom- 
pressibility coefficient ^| and is given by 



K sat (5) = 9p sat — — — \ p =p Bat . (7) 



Its value depends on the isospin asymmetry 8. Up to 4 th-o rder in 5, the isobaric 
incompressibility coefficient K sa t (S) can be expressed as ^21 

K sat (6) =K + K s&t , 2 5 2 + K s , tA 6 4 + 0{5 6 ) (8) 



with 



tf»t,a = K sym ~6L-^L, (9) 

AO 

-Ksat,4 = KsymA ~ 6-£/ S ym,4 



JoL S yniA , 9i 2 J synl L 
K K K 



IqL 2 JoK sym L 3JqL 2 JqL 2 



2K 2 K 2 K 2 2KI 



(10) 



It is interesting to see that with precision up to 4th-order in 6 for the isobaric 
incompressibility coefficient, we only need to know 8 coefficients Kq, Jo, Iq, L, 
Ksynn Jsymi ^sym,4i ^sym,4 which are defined at the normal nuclear density pg. 
Furthermore, the 4 coefficients Kq, Jo, L, and i^ S ym determine completely the 
isobaric incompressibility coefficient with precision up to 2nd-ordcr in S. 

It is generally believed that one can extract information on K s& t,2 by measuring 
the GMR in neutron-rich nuclei SBI. Usually, one can define a finite nucleus incom- 
pressibility Ka(N, Z) for a nucleus with N neutrons and Z protons (A = N + Z) 
by the GMR energy S G mr, i.e., 



Egmk — \ ^T-oT — ") (11) 

y m (r z ) 

where m is the nucleon mass and (r 2X j is the mean square mass radius of the nucleus 
in the ground state. Similar to the semi-empiric al m ass formula, the finite nucleus 
incompressibility Ka(N, Z) can be expanded as SSI 

K A {N, Z) = K + K smi A"^ 3 + K curv A~ 2 / 3 



+{ K T + K ss A-^)l^^j + 



Z 



2 



K CouiJm + --; (12) 



Neglecting the K curv term, the K ss term and other higher-order terms in Eq. (|T2|) . 
one can express the finite nucleus incompressibility K A (N, Z) as 

I ' N — Z\ 2 Z 2 
K A (N, Z) = K + K sur {A~ 1/3 + K T — + ^Coui^r, (13) 



A J ^^A^ 3 ' 
where K , K SUI {, K T , and K cou \ represent the volume, surface, symmetry, and 
Coulomb terms, respectively The K T parameter is usually thought to be equiv- 
alent to the -Ks a t.2 parameter. However, we would like to stress here that the -ftT sa t,2 
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parameter is a theoretically well-defined physical property of asymmetric nuclear 
matter as shown previously while the value of the K T parameter may depend on 
the details of the truncation scheme in Eq. (TT2J) . As shown in Ref. I^ZI, K T may 
be related to the isospin-dependent part of the surface properties of finite nuclei, 
especially the surface symmetry energy. Therefore, cautions are needed to interpret 
the K T parameter as the -Ksat.2 parameter. 

3. Results and discussions 

In the following, we show our results based on three theoretical models, namely, the 
modified finite-range Gogny effective interaction (MDI), the Skyrme-Hartree-Fock 
(SHF) approach, and a phcnomenological modified Skyrme-like (MSL) model. For 
the details of these models, we refer readers to Ref. A very useful feature of 
these models is that analytical expressions for many interesting physical quantities 
in asymmetric nuclear matter at zero temperature can be obtained, and this makes 
the study of the properties of asymmetric nuclear matter more transparent. 

3.1. Incompressibility of asymmetric nuclear matter at saturation 
density 




5 2 

Fig. 1. (Color online) The incompressibility at saturation density K SSL t as a function of <5 2 in the 
MDI interaction with x = 1, 0, and —1. Corresponding results up to 5 2 and (5 4 are also included for 
comparison. The ins ets display corresponding results at smaller isospin asymmetries with <5 2 < 0.1. 
Taken from Ref. ED 

Shown in Fig. [1] is the incompressibility at saturation density K sat (<5) as a func- 
tion of 8 2 for the MDI interaction with x = 1, 0, and — 1. It is seen that K sat {5) 
generally decreases with increasing isospin asymmetry. This feature is consistent 
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with earlier calculations based on microscopic many-body approaches Corre- 
sponding results at smaller isospin asymmetries with S 2 < 0.1 are given in the 
insets of Fig. [TJ and it shows that the result from keeping up to the S 2 term ap- 
proximates very well the exact K sa ±(S) and the contribution of if S at,4 term is not 
important as A sat (<5) displays a good linear correlation with S 2 . This feature implies 
that the absolute value of A sa t,2 is much larger than that of K aa t,i- Our results are 
also consistent with the very recent study based on the RMF model ^1 



3.2. Correlations of Jq-Kq and L-K sym 

As shown in Eq. the A sati 2 parameter is completely determined by the 4 char- 



acteristic parameters Ao, Jo, J, and A sym at the normal nuclear density. While the 
Kq parameter has been relatively well determined to be 240 ± 20 McV from the 
nuclear GMR data I " I ' ' ^ ' \ the Jo parameter is poorly known, and there is 
actually no experimental information on the Jq parameter. 
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Fig. 2. (Color online) Left panels: The Jo and the ratio Jq/Kq as functions of Kq from the 
MSL model, the MDI interaction and the SHF prediction with the 63 Skyrme interactions. Right 
panels: Correlation between K sym and L from the MSL model with 7 sym = 4/3 and 5/3, the MDI 
interaction and the SHF prediction with the 63 Skyrme interactions. The results fr om two simple 
one-parameter symmetry energies are also shown for comparison. Taken from Ref. BSl. 



Shown in the left panels of Fig. [5] are J and Jo/K as functions of A' in 
the MSL model for different values of the nucleon effective mass to* = m, 0.9to, 
0.8to and 0.7m. Also included in the figure are corresponding results from the 
MDI interaction and the SHF prediction with 63 Skyrme interactions (see Ref. 
S3 for the details). One can sec the Jq/Kq value deviates significantly from zero 
for a reasonable Ao value, which implies that the higher order Jo contribution to 
A satj 2 generally cannot be neglected. It is further seen that the correlation between 
Jo and Ao is similar among these three different models. Also, all three models 
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show an approximately linear correlation between Jq/Kq and Kq. We note that the 
correlation between Jo a nd Kq obtained in the present work is consistent with the 
early finding by Pearson While there do not exist any empirical constraints on 
the Jo parameter, we assume here the correlation between Jo and Kq from the MSL 
model is valid and then determine Jq/Kq from the experimental constraint on Kq. 

Shown in the right panels of Fig. [2] are the correlation between K sym and L 
from the MSL model with % ym = 4/3 and 5/3, the MDI interaction, the SHF 
prediction with the 63 Skyrme interactions, and two one-parameter parameteriza- 
tions on the symmetry energies. It is very interesting to see that for larger L values 
(L > 45 MeV which is consistent with the constraint from heavy-ion collision data), 
all above symmetry energy functionals from different models and paramcterizations 
give consistent predictions for the K sym -L correlation. This nice feature implies that 
using these different models and paramcterizations for the symmetry energy will 
not influence significantly the determination of the K sat ,2 parameter. On the other 
hand, the K sym -L correlation from the two one-parameter parameterizations on the 
symmetry energies deviate significantly from the MDI, MSL and SHF predictions 
for small L values. Actually, the two forms of one-parameter parametrization for 
the symmetr y en ergy may be too simple to describe a softer symmetry energy. As 
shown in Ref . 1221 the symmetry energy form E sym (p) ps 31.6(p/ J oo) 7 cannot describe 
correctly the density dependence of the symmetry energy from the MDI interaction 
with x = 1 (the Gogny interaction). Although there is no direct empirical infor- 
mation on the -Ksym parameter and some uncertainty on the K sym -L correlation 
still exist, we assume here that the correlation between K sym and L from the MSL 
model with 7 sym = 4/3 and 5/3 is valid and then use the experimental constraint 
on L to extract the value of K sym . 

3.3. Constraining the K SSLt ,2 parameter 

From above analyses, we can now extract information on the i^ S at,2 parameter from 
the experimental constraints on the Kq parameter and the L parameter. Results on 
K sa ±,2 as a function of L are shown in Figl3]for 7 sym = 4/3 (panel (a)) and 5/3 (panel 
(b)) and with K = 220, 240, and 260 MeV. In the MSL model, the correlation 
between K sym and L also depends on E sym (po) ™l To take into consideration 
of the uncertainty in the value of -E S ym(po)j we also include in Fig[3] the results 
with Kq = 220 MeV and E sym (p ) = 25 MeV as well as Kq = 260 MeV and 
E syia (po) = 35 MeV, which represent, respectively, the upper and lower bounds 
on i^sat,2 for a fixed value of L. The shaded region in Fig. [3] further shows the 
constrained L values from heavy-ion collision data, namely, 46 MeV < L < 111 
MeV. The lower limit of L = 46 MeV is obtained from the lower bound in the 
ImQMD analyses of the isospin diffusion data and the double neutron/proton ratio 
ESI while the upper limit of L = 111 MeV corresponds to the upper bound of L 
from the IBUU04 transport model analysis of the isospin diffusion 

data | 32 | 33 | 35 | 34 | 

We note that the constraint 46 MeV < L < 111 MeV obtained from heavy-ion 
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Fig. 3. (Color online) K aSy t,2 as a function of L from the MSL model with 7sym = 4/3 (a) and 5/3 
(b) and m* = 0.8m and m* = 0.7m for different values of Kq and E S ym(po)- The shaded region 
indicates the constraints within the MSL model with Kq = 240 ± 20 MeV, E S ym(po) = 30 ± 5 
MeV, and 46 MeV < L < 111 MeV from the heavy-ion collision data. The results from the SHF 
approach with 63 Skyrme interactions are al so incl uded for comparison. In addition, the constraint 
of K T = —550 ± 100 MeV obtained in Ref. p2 | 23| f rom measurement s of the isotopic dependence 
of the GMR in even-A Sn isotopes is also indicated. Taken from Ref. ESI 



collisions is consistent with the analyses of the pygmy dipole resonancesp^, the giant 

dipolc resonance (GDR) of 208 Pb analyzed with Skyrme interactions ^4zL the very 

• |49| 

precise Thomas-Fermi mode l fit to the binding energies of 1654 nuclei 1 ^, and the 
recent neutron-skin analysis^! These empirically extracted values of L represent 
the best and most stringent phenomenological constraints available so far on the 
nuclear symmetry energy at sub-saturation densities. It should be mentioned that 
the proposed experiment of parity-violating electron scattering from 208 Pb at the 
Jefferson Laboratory is expected to give an in depend ent and accurate measurement 
of its neutron skin thickness (within 0.05 fm) I 55 | 56 | anc j t impose a stringent 
constraint on the slope parameter L in future. 

From the shaded region indicated in Fig. [3l we find that for 7 sym = 4/3, we have 
-437 MeV < K Bat ,2 < -292 MeV for L = 46 MeV and -487 MeV < K satj2 < -306 
MeV for L = 111 MeV. These values are changed to -477 MeV < K sakt2 < -302 
MeV for L = 46 MeV and -461 MeV < A" sat , 2 < -251 MeV for L = 111 MeV when 
7sym = 5/3 is used. Our results thus indicate that based on the MSL model with 
4/3 < 7syra < 5/3, K = 240 ± 20 MeV, 25 MeV < E sym {p ) < 35 MeV, and 46 
MeV < L < 111 MeV, the K s&t .2 parameter can vary from -251 MeV to -487 MeV. 
The results shown in Fig. [3] is obtained from a Jq/Kq value that is evaluated with 
the default value to* = 0.8m in the MSL model. Similar analyses indicate that 
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the -Ksat,2 parameter varies from —261 MeV to —489 MeV if we use m* = 0.7m 
while it varies from —245 MeV to —485 MeV if m* = 0.9m is used. Therefore, the 
extracted value for isT sa t,2 is not sensitive to the variation of the nucleon effective 
mass. The MSL model analyses with 4/3 < 7 sym < 5/3, K = 240 ± 20 MeV, 
E sym (po) = 30 ± 5 MeV, 46 MeV < L < 111 MeV, and m* = 0.8 ± 0.1m. thus 
lead to an estimate value of i^ sa t,2 = — 370 ± 120 MeV for the 2nd-order parameter 
in the isospin asymmetry expansion of the incompressibility of asymmetric nuclear 
matter at its saturation density. 

4. Summary and conclusions 

We have studied the incompressibility of an asymmetric nuclear matter at its sat- 
uration density. We find that in its power series expansion in isospin asymmetry 
the magnitude of the 4th-order -Ksat,4 parameter is generally small compared to 
that of the 2nd-order K SSL t,2 parameter, so the latter essentially characterizes the 
isospin dependence of the incompressibility of asymmetric nuclear matter at satu- 
ration density. Furthermore, the i^ sa t,2 parameter can be uniquely determined by 
the slope L and the curvature K sym of nuclear symmetry energy at normal nuclear 
matter density and the ratio Jq/Kq of the 3rd-order and incompressibility coeffi- 
cients of symmetric nuclear matter at saturation density, and we find the higher 
order Jo contribution to -Ksat.2 generally cannot be neglected. Since there is no 
experimental information on the Jo parameter and the K sym parameter, we have 
thus used the MSL model, which can reasonably describe the general properties of 
symmetric nuclear matter and the symmetry energy predicted by both the MDI 
model and the SHF approach, to estimate the value of i^ sa t,2- Interestingly, we find 
that there exists a nicely linear correlation between K syra and L as well as between 
Jo /Kq and Kq for the three different models used here, i.e., the MDI interaction, 
the MSL interaction, and the SHF approach with 63 Skyrme interactions. These 
correlations have allowed us to extract the values of the Jo parameter and the K Bym 
parameter from the empirical information on Kq, L and £ , sym (po)- In particular, 
using the empirical constraints of Kq = 240 ± 20 MeV, E sym (p ) = 30 ± 5 MeV, 46 
MeV < L < 111 MeV and a nucleon effective mass m* = 0.8 ± 0.1m in the MSL 
model leads to an estimate of if S at,2 = —370 ± 120 MeV. 

While the estimated value of K sat .2 = — 370 ± 120 MeV in the present work has 
small overlap with the constraint of K T = — 550±100 MeV for the symmetry term in 
a semi-empirical mass fo rmula-like expansion of the incompressibility of finite nuclei 

lOOIOQI 

obtained in Rcfs. 1 1 'from recent measurements of the isotopic dependence of the 
GMR in even-A Sn isotopes, its magnitude is signi ficantly s maller than that of the 
constrained K T . Recently, there are several studies I 57 | 58 | 52 | on extracting the value 
of the -Ksa.t,2 parameter based on the idea initiated by Blaizot and collaborators 
that the values of both Kq and K satt 2 should be extracted from the same consistent 
theoretical model that successfully reproduces the experimental GMR energies of a 
variety of nuclei. These studies show that no single model (interaction) can simul- 
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tancously describe the recent measurements of the isotopic dependence of the GMR 
in even- A Sn isotopes and the GMR data of 90 Zr and 208 Pb nuclei, and this makes 
it difficult to accurately determine the value of K sat ^ from these data. Also, a very 
recent study ESI indicates that the effect due to the nuclear superfluidity may also 
affect the extraction of the -Ksat.2 parameter from the nuclear GMR. As pointed 
ou f m Ref ■ these features suggest that the K T = —550 ± 100 McV obtained 
in Ref. I ^ l ^ l may suffer from the same ambiguities already encountered in earlier 
attempts to extract the Ko and K sa t,2 of infinite nuclear matter from finite- nuclei 
extrapolations. This problem remains an open challenge, and both experimental 
and theoretical insights arc needed in future studies. 

Acknowledgments 

This work was supported in part by the NNSF of China under Grant Nos. 10575071 
and 10675082, MOE of China under project NCET-05-0392, Shanghai Rising- 
Star Program under Grant No. 06QA14024, the SRF for ROCS, SEM of China, 
the National Basic Research Program of China (973 Program) under Contract 
No. 2007CB815004, the U.S. National Science Foundation under Grant No. PHY- 
0652548, PHY-0757839 and PHY-0758115, the Welch Foundation under Grant No. 
A-1358, the Research Corporation under Award No. 7123 and the Texas Coordi- 
nating Board of Higher Education Award No. 003565-0004-2007. 

References 

1. W. Zhan et al, Int. Jour. Mod. Phys. E 15, 1941 (2006); 



2. Y. Yano, "The RIKEN RI Beam Factory Project: A status report": Nucl. Instr. 
Meth. B 261, 1009 (2007). 

3. See, e.g., |http://www.gsi.de/fair/index_e.html| 

4. See, e.g., http://ganinfo.in2p3.fr/research/developments /spiral2 

5. See, e.g., Whitepapers of the 2 007 NSAC Long Range Plan Town Meeting, Jan., 2007, 
Chicago, |http://dnp. aps.org] 

6. B.A. Li, CM. Ko, and W. Bauer, Int. Jour. Mod. Phys. E 7, 147 (1998). 

7. Isospin Physics in Heavy-Ion Collisions at Intermediate Energies, Eds. Bao-An Li 
and W. Udo Schroder (Nova Science Publishers, Inc, New York, 2001). 

8. P. Danielewicz, R. Lacey, and W.G. Lynch, Science 298, 1592 (2002). 

9. J.M. Lattimer and M. Prakash, Phys. Rep. 333, 121 (2000). 

10. J.M. Lattimer and M. Prakash, Astrophys. J. 550, 426 (2001). 

11. J.M. Lattimer and M. Prakash, Science 304, 536 (2004). 

12. V. Baran, M. Colonna, V. Greco, and M. Di Toro, Phys. Rep. 410, 335 (2005). 

13. A.W. Steiner, M. Prakash, J.M. Lattimer, and P.J. Ellis, Phys. Rep. 411, 325 (2005). 

14. L.W. Chen, CM. Ko, B.A. Li, and G.C Yong, Front. Phys. China 2, 327 (2007) 
|arXiv:0704.2340| . 

15. B.A. Li, L.W. Chen, and CM. Ko, Phys. Rep. 464, 113 (2008). 

16. D.H. Youngblood, H.L. Clark, and Y.-W. Lui, Phys. Rev. Lett. 82, 691 (1999). 

17. Y.-W. Lui, D.H. Youngblood, Y. Tokimoto, H.L. Clark, and B. John, Phys. Rev. C 
70, 014307 (2004). 




November 17, 2009 15:29 WSPC/INSTRUCTION FILE Kisolll7 



Incompressibility of asymmetric nuclear matter 11 

18. Z.Y. Ma et al., Nucl. Phys. A703, 222 (2002). 

19. D. Vretenar, T. Niksic, and P. Ring, Phys. Rev. C 68, 024310 (2003). 

20. G. Colo, N. Van Giai, J. Meyer, K. Bennaceur, and P. Bonche, Phys. Rev. C 70, 
024307 (2004). 

21. S. Shlomo, V.M. Kolomietz, and G Colo, Eur. Phys. J. A 30, 23 (2006). 

22. T. Li et al., Phys. Rev. Lett. 99, 162503 (2007). 

23. U. Garg et al., Nucl. Phys. A788, 36 (2007). 

24. G. Colo, AIP Conf. Proc. 1128, 59 (2009) |arXiv:0902.3739] . 

25. J. Aichelin and C. M. Ko, Phys. Rev. Lett. 55, 2661 (1985). 

26. C. Fuchs, Prog. Part. Nucl. Phys. 56, 1 (2006). 

27. W.D. Myers and W.J. Swiatecki, Nucl. Phys. A81, 1 (1966). 

28. K. Pomorski and J. Dudek, Phys. Rev. C 67, 044316 (2003). 

29. I. Bombaci, in Ref. ^, p. 35. 

30. S. Kubis and M. Kutschera, Nucl. Phys. A720,189 (2003) 

31. C. Xu and B.A. Li, 2009, larXiv:0910.4803l 

32. M.B. Tsang et al., Phys. Rev. Lett. 92, 062701 (2004). 

33. L.W. Chen, CM. Ko, and B.A. Li, Phys. Rev. Lett. 94, 032701 (2005). 

34. L.W. Chen, CM. Ko, and B.A. Li, Phys. Rev. C 72, 064309 (2005). 

35. B.A. Li and L.W. Chen, Phys. Rev. C 72, 064611 (2005). 

36. D. Shetty, S.J. Yennello, and G.A. Souliotis, Phys. Rev. C 75, 034602 (2007). 

37. Y. Zhang et al., Phys. Lett. B664, 145 (2008) 

38. M.B. Tsang et al., Phys. Rev. Lett. 102, 122701 (2009). 

39. W. Trautmann et al., Prog. Part. Nucl. Phys. 62, 425 (2009) 

40. Z.G. Xiao et al., Phys. Rev. Lett. 102, 062502 (2009); M. Zhang et al., Phys. Rev. 
C 80, 034616 (2009). 

41. W. Reisdorf et al., Nucl. Phys. A781, 459 (2007) 

42. W.D. Myers and W.J. Swiatecki, Nucl. Phys. A601, 141 (1996); Phys. Rev. C 57, 
3020 (1998). 

43. B.A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 

44. L. Trippa, G. Colo, and E. Vigezzi, Phys. Rev. C 77, 061304(R) (2008). 

45. A. Klimkiewicz et al. (LAND Collaboration), Phys. Rev. C 76, 051603(R) (2007). 

46. JP. Blaizot, Phys. Rep. 64, 171 (1980). 

47. JP. Blaizot and B. Grammaticos, Nucl. Phys. A355, 115 (1981). 

48. L.W. Chen, Sci. China G 52, 1494 (2009) [arXiv:0910.0086] . 

49. L.W. Chen, B.J. Cai, CM. Ko, B.A. Li, C. Shen, and J.Xu, Phys. Rev. C 80, 014322 
(2009). 

50. M. Prakash and K. S. Bedell, Phys. Rev. C 32, 1118 (1985). 

51. I. Bombaci, T.T.S. Kuo, and U. Lombardo, Phys. Rep. 242, 165 (1994). 

52. J. Piekarewicz and M. Centelles, Phys. Rev. C 79, 054311 (2009) [arXiv:0812.4499] . 

53. J.M. Pearson, Phys. Lett. B271, 12 (1991). 

54. M. Centelles, X. Roca-Maza, X. Vinas, and M. Warda, Phys. Rev. Lett. 102, 122502 
(2009). 

55. C.J. Horowitz, S.J. Pollock, P.A. Souder, and R. Michaels, Phys. Rev. C 63, 025501 
(2001). 

56. R. Michaels, P.A. Souder, and CM. Urciuoli, 2005, URL 
http: / /hallaweb . j l abTorg/parity /prex | 

57. H. Sagawa, S. Yoshida, CM. Zeng, J.Z. Gu, and X.Z. Zhang, Phys. Rev. C 76, 
034327 (2007). 

58. J. Piekraewicz, Phys. Rev. C 76, 031301(R) (2007). 

59. E. Khan, Phys. Rev. C 80, 011307(R) (2009). 



